***********************************************************************
***********************    HETEROGENEITY    ***************************
***********************************************************************
scalar list 
gen P_uvalue256_fem = _bAR_fem_a*a + _bAR_fem_b*b +  _bAR_fem_c*c + _bAR_fem_d*d + _bAR_fem_e*e + _bAR_fem_f*f +  _bAR_fem_g*g + _bAR_fem_h*h
gen I_uvalue256_fem = _bHE_fem_a*a + _bHE_fem_b*b +  _bHE_fem_c*c + _bHE_fem_d*d + _bHE_fem_e*e + _bHE_fem_f*f +  _bHE_fem_g*g + _bHE_fem_h*h

gen P_uvalue256_men = _bAR_men_a*a + _bAR_men_b*b +  _bAR_men_c*c + _bAR_men_d*d + _bAR_men_e*e + _bAR_men_f*f +  _bAR_men_g*g + _bAR_men_h*h
gen I_uvalue256_men = _bHE_men_a*a + _bHE_men_b*b +  _bHE_men_c*c + _bHE_men_d*d + _bHE_men_e*e + _bHE_men_f*f +  _bHE_men_g*g + _bHE_men_h*h

gen markerM = "♂"   // Male symbol
gen markerF = "♀"   // Female symbol

twoway ///
	(scatter P_uvalue256_fem I_uvalue256_fem ,   msize(tiny) mcolor(orange%60) mlabel(markerF) mlabcolor(orange%60)  mlabposition(0)) ///
	(scatter P_uvalue256_men I_uvalue256_men ,   msize(tiny) mcolor(ebblue%80) mlabel(markerM) mlabcolor(ebblue%80) mlabposition(0)) ///
	(scatter P_uvalue256_men I_uvalue256_men if cardname=="00000000",   msize(tiny) mcolor(black%80) mlabposition(0)) ///
, ///
xline(0, lpattern(-) lc(gray)) ///
yline(0, lpattern(-) lc(gray)) ///
yt("Palestinians `Acceptability'") ///
xt("Israelis `Acceptability'")  ///
text(-0.04 0.1 "status quo", size(vsmall)) ///
note(" " "", size(vsmall)) ///
legend(off)  ///
graphregion(fcolor(white)) 
graph display, ysize(1) xsize(1)